Can Non-Linear Elasticity Explain Contact-Line Roughness at Depinning? 
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We examine whether cubic non-linearities, allowed by symmetry in the elastic energy of a contact line, may 
result in a different universality class at depinning. Standard linear elasticity predicts a roughness exponent 
£ = 1/3 (one loop), £ = 0.388 ± 0.002 (numerics) while experiments give £ ~ 0.5. Within functional RG we 
find that a non-local KPZ-type term is generated at depinning and grows under coarse graining. A fixed point 
with £ w 0.45 (one loop) is identified, showing that large enough cubic terms increase the roughness. This fixed 
point is unstable, revealing a rough strong-coupling phase. Experimental study of contact angles 9 near n/2, 
where cubic terms in the energy vanish, is suggested. 



Experiments measuring the roughness of the contact line of 
a fluid wetting a disordered solid substrate have consistently 
found a value ( « 0.5 for the roughness exponent. This re- 
sult is highly reproducible from superfluid Helium 1 1] to vis- 
cous glycerol-water mixtures |2], and in situations which can 
rather convincingly be argued to be at or at least very near the 
depinning transition. Explaining this high value for ( poses a 
theoretical challenge. It may result in a broader understand- 
ing of the depinning transition in other systems, since similar 
values are also measured in cracks |3]. 

The simplest elastic model of a contact line |4] consists of 
an effective elastic energy, quadratic in the height field h(x) 
(displacement in the solid plane), with non-local dispersion 
c\q\ (long-range elasticity) due to the surface tension of the 
fluid meniscus. Substrate inhomogeneities are modeled by a 
random-field disorder coupling to h(x). The resulting model 
for the depinning transition of an elastic manifold (general- 
ized to d internal dimensions, here d = 1) has been exten- 
sively studied, and the predictions debated for some time. 
Functional RG (FRG) methods were developed initially to one 
loop 1 5] predicting ( = e/3 to all orders, here e = 2 — d, 
identical to the statics of random field. Careful analysis be- 
yond one loop however revealed new irreversible terms in the 
RG which clearly distinguish statics and depinning, and yield 
C = e/3(l+0.397e)+O(e 3 ) Jj|. Novel high-precision numer- 
ical algorithms found 0| < = 0.388 ± 0.002 midway between 
the one- and two-loop results. This value is too low to account 
for the experiments. 

Various mechanisms have been proposed |8] such as lat- 
eral waves or plastic-type dynamics. It is unclear whether any 
of them are universal enough to explain the robustness of the 
experimental values for £. More complex dissipation mecha- 
nism may be at play but they should not be important for the 
roughness if, as believed |2], the experiment is at quasi-static 
depinning. 

Before abandoning the elastic model, one must first check 
for neglected effects. It has been known for some time, for 
conventional linear elasticity cq 2 , that there is another uni- 
versality class, anisotropic depinning SHU [Till . There, a 
non-linear Kardar-Parisi-Zhang (KPZ) term A(Vft-) 2 becomes 
relevant, resulting in a singular dependence of the threshold 
force on rotation of the contact line in the plane. That such 
a dissipative term is generated in the equation of motion at 



velocity v > is straightforward, and it was shown recently 
within the FRG that it survives 11211 even at quasi-static de- 
pinning v — > + , only if anisotropy exists in the substrate 
disorder or in the motion of the manifold. Although Ref. fLlll 
shows that a larger exponent is possible for long-range elas- 
ticity, an explanation based on this term alone is problematic. 

An important feature of the contact-line problem is that at 
contact angles different from 8 = ir/2 the symmetry h — > —h 
is absent. Thus the elastic energy contains cubic non-linear 
and non-local terms. The equation of motion thus contains 
non-linear terms breaking this symmetry, even in the absence 
of a driving force fl3ll . Being different in nature from the con- 
ventional dissipative KPZ term, it is important to understand 
their effect at quasi-static depinning. Their effect in the mov- 
ing phase has been recently investigated and found to 
lead to a roughening transition analyzed in connection with 
the development of a Landau-Levich film. 

The aim of this letter is to examine whether non-linear 
terms in the energy may result in a different universality class 
at depinning. Although naively irrelevant, they do generate 
at depinning a non-local KPZ-type term which grows under 
coarse-graining. We find via a FRG calculation two possi- 
ble phases separated by a fixed point at a critical value of the 
disorder. Interestingly, the roughness exponent at this critical 
point is, within one-loop accuracy, ( w 0.45 > 1/3. Al- 
though we do not control the rough phase at strong disorder, 
this shows that ( can be increased by non-linear terms. This 
growing non-local term also arises in the moving phase, with 
similar features and a significantly smaller exponent. We dis- 
cuss the interest of studying the system for interfaces with 
contact angles in the vicinity of 8 = tt/2, the point where 
cubic terms in the energy vanish. 

Let us describe the model of Ref. 0,[3 for a fluid wetting 
a flat solid in coordinates suitable for any (global) equilibrium 
contact angle 8 = 8 e , see Fig.^ The liquid-air interface (LA) 
is denoted by {x, y, z(x, y)} (flat when z(x, y) = 0), the flat 
solid (S) surface by {x, y, —y tan 8}. They meet at the contact 
line y = h(x) cos 8, included in (S), the boundary condition 
z(x,y = h{x) cos 8) = —h(x) sin 8. The total energy is E = 
E LA + E SL with: 



Ela [z] = 7la / da; 



y>h(x) cos 9 



dyA[z(x,y)] (1) 
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FIG. 1: The geometry of the contact-line. 



£sl= fdxf dy'j(x,y'), (2) 

J Jy'>h(x) 

where A = \Jl + (d x z) 2 + (d y z) 2 gives the LA area, and the 
wetting area energy density 7 = 7sl — 7as is a random func- 
tion of the in-plane position (y' = yj cos 8). Minimizing the 



LA interface energy E-^a for fixed h(x) yields the equilibrium 
profile y), and one can show that the relation between 
force and contact angle holds locally, SEi J A[zh]/Sh(x) = 
7la cos 6(x). One defines the equilibrium contact-angle value 
6 e via the average force, i.e. through 7(2;, y 1 ) — 7la cos# e + 
j(x,y') with 7(1,?/) = 0. 

Expansion of £la[z/i] up to third order in h proceeds by 
solving V 2 Zh = in the form Zh(x,y) = f a q e lqx ^^ v . 
Solving for the boundary conditions and inserting in E^a [zh] 
yields, up to terms of 0(h 4 ): 

E hA [h] = ^ j '\q\h q h- q + - J [\q\\k\ + qk]h k h q h_ q _ k 

(3) 

with ci = 7sin 2 (0), A = c\ cos{6). The form of the cu- 
bic term e quq2tq3 h qi h q2 h q3 , J2<H = °> witn e ?i, «2,? 3 ~ 
^2i<j i s tne on iy possibility imposing that e is 

symmetric, homogeneous of degree 2, and vanishes for qi = 
(invariance under a uniform shift h(x) — > h(x) + est). As ex- 
pected the symmetry h k — > — is restored for 8 = ir/2 
where A = 0. To the same order the general form of the equa- 
tion of motion (EOM) is: 



r]9thk,t = -c(k)hk,t ~ 2 J ^ 2 ^' k ~ ^ + ^ 3 (~ fc ' ?) + 



where the pinning force has correlator F{x, h)F{x' , h') = 
S d (x — x')A(h — h'), and a thermal noise may be added. 
In view of later RG calculations, we slightly generalize the 
model, defining (with normalized vectors q = q/\q\): 

*i(qi,<h) = \qi\ a \<h\ a fitfi-&), c(q) = c a \q\ a (5) 

The contact line corresponds to a = 1. For reasons detailed 
below we consider the parameterization: 

f 2 (z) = X 3 (l + g(z)) + X 0l f 3 (z) = X 3 (l + g(z)) (6) 

where g(z) — z or more generally an odd function with 
g{—l) = — 1. This EOM is obtained from in the sim- 
plest case, which assumes fast relaxation of the meniscus and 
dissipation via molecular jumps Ed: 

7}dth{x,t) -8E[h) 

y/l + ld x hY ~ Sh(x,t) 

= 7[cos 8(x, t) — cos 6 e ] + j(x, h(x, t)) , (7) 

where 77 is a dissipative coefficient (note that the above equa- 
tion neglects viscous hydrodynamic losses inside the moving 
liquid wedge | 

At zero or vanishingly small velocity, i.e. at the depinning 
threshold, the non-linearity (d x h) 2 on the left-hand side can 
be neglected to this order, and using the form Q one finds 
and © with A 3 = A, A = and F(x, h) = 7(2;, h). 
These are the microscopic (bare) values, and we show below 
that the form given in (|3J and is preserved under RG at 



k,k-q)]h q , t h k - q>t + / [F(x, h{x) + vt) + f - t]v] e lhx 

J X 

(4) 

I 

depinning. The crucial point is that, while forbidden at the 
bare level by the potential form of the EOM rjdfh = —5E/8H, 
a non-zero and positive value for Ao will be generated beyond 
the Larkin length from the non-analyticity of the renormalized 
force correlator. It corresponds to the generation of a non- 
local KPZ-type term. The case of more complicated dynamics 
is mentioned below. We now focus on the analysis of @. 

FRG starts by calculating the 1 -loop corrections to (0} and 
its associated dynamical action, perturbatively in A(h) and 
the Xi. Graphical rules are given in Fig. [2] The graphs rep- 
resenting the corrections to c Q , Xi, A(h) and rj are shown in 
Fig.0 To illustrate the main ideas, we explain how the first 
correction to c a is computed. It is the sum of two contribu- 
tions (see Fig. 0J: 

-A'(0+) f \k\ a \ P \ a r 

(8) 

where q = —(p + k). These are depicted on Fig. |4]in the order 
of their appearance in (|8j. Denominators originate from time 
integrals of the bare response R kt = (rj)~ 1 e~ Ca ^ k ^ t ^ 1 Q(t), 

~ 1 1 >* 

FIG. 2: The graphical rules: Propagator, disorder vertex, non- 
linearity ~ A2 and non-linearity ~ A3. 
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FIG. 3: The diagrams correcting c a (first line), A; (second line), 
A(/i) (third line) and r\ (last line). Only diagrams proportional to 
Ai are shown. 



while numerators come from the cubic vertices. The factor 
A'(0+) is A'(v(t - t' j), taken in the limit of v -> 0+ and is 
non-vanishing only when the cusp is formed, i.e. beyond the 
Larkin length. Expanding in vanishing external momentum p 
yields: 



Sc 



A'(0 H 



1 



(Mz) + M-z)) 



(9) 



where (f(z)) denotes the angular average ~ J_. dz (1 — 

z 2 )^' j '(z) (with (1) = 1). Appearance of the combination 
(f2(z) + fa{—z)) — 2A3 + Ao for arbitrary odd g(z) in all 
graphs is a general feature and yields the (re)definition of the 
vertex explained in Fig.|4] It gives the "anomalous" correction 
(i.e. resulting from the cusp) to c a to which the second one of 
Fig.|5]must be added. One shows that there is no correction 
to the non-linear vertex function fs(z) because that would in- 
volve necessarily another vertex with an external h leg at 
zero external momentum, yielding fs(— 1) = (which re- 
mains true since there is no mechanism to correct it). Thus 
one finds <5A3 = 0. In the graphs correcting /2 (second line of 
Fig. |3} appearance of the above-mentioned combination im- 
plies that only the uniform part of j^, i.e. Aq is corrected: 



<5A - (2A 3 + Ao) 2 



A'^ 



3— i^-Ao 



(10) 

While in standard (thermal) KPZ the two corresponding di- 
agrams have opposite sign because of the derivatives and A 
is uncorrected (a consequence of Galilean invariance), here 
the presence of absolute values of the momenta results in the 
same sign. 

Corrections to disorder (line 3 of Fig. [5J are the same as 
in O- Power counting then leads to the definitions A(h) = 



Af- 2a+2C c~ 



A(hAp), A, = Af^c^Xi, where A/ = A 



is the running UV cutoff (e.g. in a Wilson approach). One 
finally obtains the following set of FRG equations (with e = 





FIG. 4: The first diagram from Fig.[3]correcting c a , see |8j. 



2a - d): 

di lnc Q = A'(0 + )(2A 3 + A Q ) - A(0)(2A 3 + A )A 
di In A 3 = C - a - di In c a 

9 ; A = (C - a - di lnc Q )A - 2A'(0+)(2A 3 + A ) 2 

+3A(0)Ao(2A 3 + A ) 2 (11) 
diA(h) = (e - 2( - 2di hxc a )A(h) + (hA'(h) 

+\~A{h) 2 Xl - lA'(h) 2 + A"(h)(A(h) A(0))] 
a i ln7 ? =iA'(0+)A -A"(0+). 

They admit the standard attractive depinning fixed point cor- 
responding to linear elasticity (isotropic IS depinning class): 
A3 = Ao = and Aj S (/i) with £ = e/3 to this order (cor- 
rected at two loops). To order e, this FP is stable to adding 
a small non-zero A3 , A which have linear eigenvalue ( — a. 
This FP controls a phase with small Ao and A3. Ao is gen- 
erated from A3 beyond the Larkin length and the ratio A0/A3 
goes to a constant in this phase. 

We found a second fixed point which controls at given dis- 
order the transition between the small Ao phase and a large Ao 
regime (strong coupling). One easily sees that the ratio A3/A0 
flows to zero at this transition. To look for the FP one can thus 
set A3 = and redefine A(h) — AqA(/i/Aq), yielding: 



diA(h) = [e - 2C - 2A'(0+) + 
+ l -A{h) 2 [A'{hf 



2A(0)]A(h) + (hA'(h) 

+ A"(h)(A(h)-A(0))} 
(12) 



with di lnAo = C - C and C = a + 3A'(0+) - 4A(0). The 
FP function obtained numerically A*(h) = ef a (h) is posi- 
tive and short-ranged, and f a depends only on a = a/e; the 
associated roughness exponent is 



C = 0.450512e 



(13) 



for a = 1, which should be compared to the value (is = e/3 
to the same 1-loop accuracy, demonstrating the increase in the 
roughness exponent due to non-linearities. If we assume the 
same relative increase of ( due to the higher-loop corrections 
0] from 1/3 to the observed Cis ~ 0.388 ± 0.002 0], we 
would obtain here £ « 0.53, tantalizingly close to the exper- 
iments. The dynamical exponent is z — a + di \n(rj/c), i.e. 
z = 1 - 0.205213e yielding z = 0.79487 for the physical 
case. As for anisotropic depinning, a third exponent is nec- 
essary here, ip — — dilnc a — 0.170449e, and scaling yields 
the correlation length exponent v = l/(a — C + an d tne 
velocity-force exponent (3 = v(z — £) w 0.478 such that 
v ~ (f ~ fc) 13 - This fixed point is unstable in one direction 
(leading eigenvalues /ii = 0.938 and //2 = —1.23) consistent 
with the existence of two phases. The strong-coupling phase 
cannot be accessed by the present method, but the increase of 
C is likely to persist there. Since this FP is attractive in all 
other directions the experiments may be susceptible to a very 
long crossover dominated by this FP. This can be tested by 
careful numerical integration of di Q , not attempted here. 

The FRG equations di It have been derived within a 
double expansion in e and a. The question arises of 
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whether operators with more non-local derivatives |V| Q /i ~ 
\k\ a h are indeed irrelevant, as suggested by power-counting. 
A detailed analysis shows that in the space of perturba- 
tions where one adds to A(h x t — h x y) the two func- 
tions A 2s (h xt - h xV ) (\V\ a h xt + M a h xt >) + A 2u (h xt - 
h X f) (\V\ a h x t — \V\ a h X f), the largest eigenvalue is —0.12 
(for e = a = 1), indicating no additional instability of the 
FP. These terms arise at the bare level for a more complicated 
dynamics, but should not change the result in the quasi-static 
limit studied here 1 14[ • 

We now sketch the analysis of the moving case, very near 
depinning, v > small, for details see 1 14]. Since the gen- 
eration of Ao is a new feature, we reexamine [13]. At large 
scales in the moving phase the quenched pinning force acts as 
a (thermal) white noise of strength 2D (notation as in lfl311 ) 
and thus r\ is uncorrected. Using the same parameterization as 
above we find H13I1 : 



(14) 

2 



di lnc Q = z — a — tS^I + 2r) 

di In A = z + ( - 2a + -g(l + 2rf 



dig=-{d + a)g + g 2 



\ + \{l + 2r)+ l -(l + 2r) 

with dilnr — — \g{l + 2r) 2 and we have defined c 

c a A l , A, = AjA ; 
g = c~ 3 D\q, and r 



2C.-Z 

A3/A0. These equations exhibit a 



weak-coupling phase controlled by the g = attractive fixed 
point (which, in the physical case corresponds to logarithmic 
roughness Q — 0) and a strong-coupling phase. They are sep- 
arated by a FP at g = g* = jy (d + a), with universal values 
for the exponents z — a + (d + a) and £ = ^ 
i.e. for a = d = 1: 



3 d+a 
11 2 



C « 0.273 



1.364 



(15) 



This rather low value for £ suggests that a scenario based on 
a slowly moving contact line is not adequate to explain the 
experimentally observed roughness. 

In conclusion, we have examined using RG the effect of 
non-linear elasticity for the contact line at and near depin- 
ning. We found that even for isotropic disorder, a non-local 
KPZ term is generated and may, for large enough bare non- 
linear elasticity and disorder, destabilize the standard linear- 
elasticity depinning fixed point, yielding values for the rough- 
ness exponent compatible with experiments. This scenario 
could be tested by high-precision numerics. It may also be 
explored experimentally by carefully choosing the fluid and 
the solid substrate fl8il : since all odd non-linear terms in the 
elastic energy vanish at 9 = ir/2 one can surmise that the total 
effect of non-linear terms, and hence the apparent contact-line 
roughness, is minimal there QsJ. 

We thank C. Guthmann, W. Krauth, S. Moulinet, E. Rolley, 
A. Rosso and T. Vilmin for interesting discussions. 
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